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C/^ ' Abstract 



We derive a modified Buchdahl inequafity for scalar-tensor theories of gravity. In gen- 
eral relativity, Buchdahl has shown that the maximum value of the mass-to-size ratio, 
2M/R, is 8/9 for static and spherically symmetric stars under some physically reasonable 
assumptions. We formally apply Buchdahl's method to scalar-tensor theories and obtain 
theory-independent inequalities. After discussing the mass definition in scalar-tensor the- 
' ories, these inequalities are related to a theory-dependent maximum mass-to-size ratio. 

, We show that its value can exceed not only Buchdahl's limit, 8/9, but also unity, which 

' we call the black hole limit, in contrast to general relativity. Next, we numerically examine 

' the validity of the assumptions made in deriving the inequalities and the applicability of 

00 , our analytic results. We find that the assumptions are mostly satisfied and that the mass- 

0^ ' to-size ratio exceeds both Buchdahl's limit and the black hole limit. However, we also 

"'^ I find that this ratio never exceeds Buchdahl's limit when we impose the further condition, 

p — 3p > 0, on the density, p, and pressure, p, of the matter. 



1 Introduction 



Einstein's general relativity postulates that gravitational interactions are mediated by a 
tensor field, g^^^. It is also well-known that electro-magnetic interactions are mediated by a 
vector field, A^. One may therefore suspect that some unknown interactions may be mediated 
by scalar fields. Such theories have been suggested since before the appearance of general 
relativity. Moreover, it has been repeatedly pointed out over the years that unified theories 
that contain gravity as well as other interactions naturally give rise to scalar fields coupled 
to matter with gravitational strength. This motivation has led many theoretical physicists to 
study scalar-tensor theories of gravity (scalar-tensor theories) The scalar-tensor 

theories are natural alternatives to general relativity, and gravity is mediated not only by a 
tensor field but also by a scalar field in these theories. Recently, such theories have been of 
interest as effective theories of string theory at low energy scales 

Many predictions of the scalar-tensor theories in strong gravitational fields are summarized 
in Ref.Q,p,[Q. It has been found that a wide class of scalar-tensor theories can pass all 
experimental tests in weak gravitational fields. However, it has also been found that scalar- 
tensor theories exhibit different aspects of gravity in strong gravitational fields in contrast to 
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general relativity. It has been shown numerically that nonperturbative effects in the scalar- 
tensor theories increase the maximum mass of an isolated system such as a neutron star[^,[Q. 

In general relativity, the mass-to-size ratio of a star has physical significance, especially for 
an isolated system. Buchdahl has obtained a maximum value of the mass-to-size ratio of a 
static and spherically symmetric star under the following physically reasonable assumptions 




• No black hole exists. 

• The constitution of the star is a perfect fluid. 

• The density at any point in the star is a positive and monotonously decreasing function 
of the radius. 

• An interior solution of the star smoothly matches an exterior solution, i.e., Schwarzschild's 
solution. 

Buchdahl has obtained an upper limit of the mass-to-size ratio as 2A//i? < 8/9. We shall refer 
to this as the Buchdahl inequality. 

Motivated by Buchdahl's theorem, we shall derive a modified Buchdahl inequality to obtain 
the maximum mass-to-size ratio in scalar-tensor theories. We then numerically examine the 
validity of the assumptions made in deriving the inequality. The applicability of our analytic 
results is also examined. This paper is organized as follows. In section 2, we summarize the 
basic equations in the scalar-tensor theories. In section 3, we derive a modified Buchdahl in- 
equality in the scalar-tensor theories, and the numerical results are compared with the analytic 
results in section 4. A brief summary is given in section 5. 



2 Basic equations 



We shall consider the simplest scalar-tensor theory [0,0, flltl . In this theory, gravitational 
interactions are mediated by a tensor field, g^i^, and a scalar field, 0. The action of the theory 
is 



S = 



1 

16^ 



0^ 1—9 9,ti(P,v 



d X -\- Sraatter\^ rm 9nij\^ 



(2.1) 



where ^(0) is a dimensionless arbitrary function of 0, represents matter fields, and Smatter 
is the action of the matter fields. The scalar field, 0, plays the role of an effective gravitational 
constant as G ~ 1/0. Varying the action by the tensor field, g^jy, and the scalar field, 0, yields, 
respectively, the following field equations: 



Guv 



Stt 



c.(0) 



1 



3 + 2w(0) 

Now we perform the conformal transformation, 



(^^,^.<P.. - ^5p.3"''0,a0,/3) + ^(V^V,0 - g^.n0), (2.2) 

(2.3) 



8^r - ^5"'50,„0.0 



(2.4) 
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such that 



G.A^i^) = \, (2.5) 



where G, is a bare gravitational constant, and we call A{(p) a coupling function. Hereafter, 
the symbol, denotes quantities or derivatives associated with g*^^. Then the action can be 
rewritten as 

S = — ^ / V^*{R* - 2g^''ip^^^^^)d*X + SmaUer[-9m,A\ip)g,^,], (2.6) 



167rG 

where the scalar field, ip, is defined by 

Varying the action by g^^^i, and ip yields, respectively, 

G^fji, = SttG^T*^^ + 2 (^(p,fj.ip,i, - ^g*i_t,,g^'^ip,a'p,p^ , (2.8) 
D^(p = -47rG,a(^)T,, (2.9) 
where Ti*" represents the energy-momentum tensor with respect to defined by 

jifiv _ 2 SSmatter [^m ; ('i'^).9*pi^] _ 74^((p)T'^'' (2 10) 

V-9* Sg^^^i, 

The conservation law for T^'^ is given by 

y*uT:^ = a((/j)T,V,^(/p. (2.11) 

The field equations (2.9) and ( ^.ll| ) tell us that the coupling strength, oi{if>), plays a role 
in mediating interactions between the scalar field, if, and the matter. General relativity is 
characterized by having a vanishing coupling strength: (y^{ip) — 0, i.e., A{ip) ~ 1. The Jordan- 
Fierz-Brans-Dicke theory is characterized by having a (/s-independent couping strength: a{ip) — 
ao = const., i.e., A{ip) = e""*^ 0,1^ • Observational constraints on the couphng strength are 
summarized in Appendix A. 



3 Modified Buchdahl's theorem in scalar-tensor theories 



In this section, we consider a static and spherically symmetric space-time with a perfect 
fluid. First, we derive a modified Buchdahl inequality. Then the inequality is reformulated to 
obtain the maximum value of the mass-to-size ratio in the scalar-tensor theories. Hereafter, 
we refer to (g^v, 4>) and (5*^1/, if), respectively, as the physical frame and the Einstein frame. 
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3.1 A modified Buchdahl inequality in scalar-tensor theories 



In the Einstein frame, a line element of the static and spherically symmetric space-time is 
written as [|lO| 

dsl = -f,{r)dt^ + K{r)dr^ + r^dQ^. (3.1) 

The stress-energy tensor for the perfect fluid in the Einstein frame is given by 

Tr = [p* + P*)« + p,g^^ , u,o. = -^/U^ {dt)o,, (3.2a) 

where is the four velocity of the matter. The stress-energy tensor for the perfect fluid in 
the physical frame is given by 



(3.2b) 



where u" is the four velocity of the matter. Fluid variables in the physical and Einstein frames 
are related according to 



P = P*A^'^{'P), 
p = p^A^^{(p). 



Now the field equations (p.q) and (2.9) are reduced to the following equations: 



(3.3) 
(3.4) 
(3.5) 



r 2/, 2K 2/* r 2/i, 



= SnGt^pt^hf: + ip 



'2 



B7rG*p*/i, — (fi 



'2 



(3.6) 
(3.7) 

(3.8) 



= 47rG,a((^)(p, -3p,), (3.9) 



where a prime denotes differentiation with respect to r. As is often done in the cases of general 
relativity, we define a mass function, m^,{r), in the Einstein frame as follows: 



Then 



where 



is rewritten as 



Kir) 



1 - 



2m* (r) 



Pa!! = P* + 



'2 



That is, Pctf plays the role of an effective density in the Einstein frame. 



(3.10) 



(3.11) 



(3.12) 
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In order to derive a modified Buchdahl inequality, we assume 



K{r) 


> 


0, 


Mr) 


> 


0, 




> 


0, 




< 


0. 



(3.13) 
(3.14) 
(3.15) 
(3.16) 

These assumptions imply the following: 

• No black hole exists in the Einstein frame. 

• The effective density, p„tt, is a positive and monotonously decreasing function of the 
radius. 

Moreover, we assume that an interior solution of the above field equations smoothly matches 
the corresponding exterior one. Note that these assumptions are concerned with the unphysical 
variables and that their validity should be examined. This will be done later. 
Using the assumption ( 3.16| ), it is easy to verify the following inequality: 



(?y^»- 



Moreover, with (|^, (|U) and (|1|), we obtain 



(3.17) 



\/y*h* 




> 0. 



(3.18) 



Accordingly, we have 



ri < r-i. 



(3.19) 



Now let the inequality ( 3.19| ) be reformulated in terms of variables of the exterior solution. 
The exterior solution, whose derivation is given in Appendix B, is 



where 



eA(x) 




ax, 


e7(x) 




--) 

xj 






^In 

a 



X 



(3.20) 

(3.21) 
(3.22) 

(3.23) 



Here a, 5, c and fo are constants of integration, and (po is the asymptotic value of if at infinity. 
Moreover, the constants, a, b and c, must satisfy the following relation (Appendix B): 



4c" 



(3.24) 



One may expect that % = a is an event horizon. However, this is not the case in generic 
scalar-tensor theories, where the null surface, x = a, is a curvature singularity in the Einstein 
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frame. The singular nature of the unphysical space-time at x = a can also be seen when 
transformation to the Schwarzschild coordinate is made. The Schwarzschild coordinate, r, and 
the Just coordinate, x, are related by the following relation: 



r X 1 



(3.25) 



Xs 1 - — 

Xs 



Mrs) = 1- — 

Xs 



K{r,) = 1 1 - 



a + b 



One finds that, when a 7^ 6, x — ^ the Just coordinate corresponds to r = in the 
Schwarzschild coordinate. 

By matching the interior solution to the exterior solution, we obtain the following relations: 



(3.26) 
(3.27) 
(3.28) 

(3.29) 
(3.30) 
(3.31) 



where the subscript, s, refers to x evaluated at the surface, r = r^. Note that 

Xs> a> b. 

Since h^{r)~^ = 1 — 2m^{r)/r, ( 3.28| ) becomes 

(o + b) 



4X. 



a 

1 

Xs 



Accordingly, by virtue of the positivity of (r) , we obtain the additional inequality 

4x. 



With ( 3.19 ) and (3.25) ^ (3.28), we obtain the following relation for r < 



Integrating (3.32) from the center, r = 0, to the surface, r — Vs, we obtain 



< Vim 



^' s Jo 



< Vf*{rs) 



Or 3 , 



r\/KAr)dr 

2m, (rs) 



r 1 



dr 



2r7i,(rs) 



4m, (r^) 



(3.32) 



(3.33) 



where the inequality (3.17) has been used. 
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The inequalities obtained to this point can be simpHfied in terms of new parameters defined 



by 



a h 
— , bs = — , 

Xs Xs 



Substituting ( ^!3C| ) into (|!3|), we obtain the foUowing inequahty: 



< V^(o) 



< -^{l-as 



{as + hsY 



2hs-^^^^±^^hsVl^s 



I- ^ I 



46. 



{as+bsf 



4(1 -a,) 



(3.34) 



(3.35) 



The above inequahty can be further simphfied, and we finahy obtain a modified Buchdahl 
inequahty in the scalar-tensor theories as 



F{as,bs 



36, 



1 



{as + bs){as + 26s) - bs\/\ - as > 0, 



(3.36) 



supplemented with (3.29) and (3.31). 

The modified Buchdahl inequality can be solved to yield 



bs < as < 2^67 - bs 
bs<as < 2V26I- 26s 
Forbidden 



for < 6s < 4(3- 2\/2), 
for 4(3-2^2) < 6s < |, 
for 6s > I. 



(3.37) 



Fig. 1^ displays the allowed region, D , of ( qs, 6s). 

In addition, we can obtain, with ( 3.24 ), the following upper limit on |cs 



Cs < 



(3.38) 



The inequality ( p!37| ) is significant. The third inequality of ( |3.37| ) gives us a necessary 
condition for a spherical star to exist and is reduced to Buchdahl's theorem in general relativity 
when c = 0, i.e., a = b, and, accordingly, x = r. In this case, we have (i? = Ts) 







2M 8 
~R~ - 9' 



(3.39) 



where M is the ADM mass defined at spatial infinity, and R is related to the surface area, S, 
as 5 = AttR^. 

The new and important inequality ( 3.38 ) is characteristic of the scalar-tensor theories and 
does not have a general relativistic counterpart. It has been found the appearance of non- 
perturbative behavior of the scalar field in the previous numerical studies [11,0. Our result 
implies that, even in a strong gravitational field, the characteristic amplitude of the scalar field, 
\cs\, is bounded. 

It is important to note that we have not used any assumption regarding the coupling 
function, A(ip), in deriving the inequalities. In particular, the inequalities (3.37) and ( 3.3^ 
give theory- independent constraints on the parameters, 6s and Cs- 
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3.2 The mass-to-size ratio 



Now we reformulate the inequalities derived in the previous section, which are in terms of 
variables in the Einstein frame, in order to obtain the mass-to-size ratio in the physical frame. 
To do this, the coupling function, A(ip), should be specified. In this paper, we assume as an 
example of this coupling function the simple form 



A{ip) = e^^^' 



where /9 is a constant |^ , . Then the coupling strength, a{(p), becomes 

a{ip) = (3(p. 



(3.40) 



(3.41) 



A natural definition of the radius of a spherical star is obtained by using its (physical) 
surface area as follows. In the physical frame, the surface area, S, is given by 



47rv42(<^^)e^(x3)-7(xs) 



exp 



/3 liln(l_a,) 



(3.42) 



where we take the asymptotic value of the scalar field as ipo = 0, and, accordingly, we have 
A{ipo) = 1 and a{ipo)=0. This surface area defines the physical radius, R, of the star in a 
similar manner as in general relativity: 



V 47r 



(3.43) 



When ipQ — 0, the effective gravitational constant, G, defined in Appendix A, is equal to G,. If 
ipo ^ 0, contributions of the scalar field appear in the above expression of R in terms of c (3 ipo. 

The definition of the mass in the Brans-Dicke theory is found in Ref.[|l^, and the mass 
in the scalar-tensor theories is defined in the same manner. In defining the mass, the metric 
should be expressed in the isotropic coordinate, f, which is related to x 



In our model, the exterior solution is then rewritten in terms of f as 
'1 - - 



(3.44) 



1 + — 



dr 



G^(j) = A ^(i/?) = exp 



4/5c2 



1 - — 
In— ^ 

1 + 



(3.45) 



By introducing the asymptotic Cartesian coordinates such that r = a/ {x'^Y + {x^Y' + (a^^)^) 
the asymptoric form of the solution is easily found to be 



G*( 
9^j 



0.1 

r 



2Ms 



900 



-1 + - = -1 
r 



2{Mt + Ms) 



6^J = I 1 



2(Mt - Ms) 



(3.46) 
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where the quantities, Ms and Mt, are called, respectively, the scalar mass and the tensor mass. 
At Newtonian order, their sum, M = Mt + Ms, plays the role of the mass and is called the 
active gravitational mass. In our model, Ms = 0, and, accordingly, Mt = M = 6/2. Hereafter, 
we call M the mass for simplicity. 

Now we are ready to calculate the mass-to-size ratio in the scalar-tensor theory as a function 
of Os, bs, Xs and a specific parameter of our model, /3. We obtain 

b bs-as 

H{as, bs; /3) = — = 6^(1 - a^) exp 
K 

In Fig.^, in the allowed region of (a^, bs), we display lines on which H(as, bs] /?) is equal to 8/9 
for various values of /?. For a fixed value of f3, the region above the line corresponds to the case 
that the mass-to-size ratio, 2M/R, exceeds Buchdahl's limit, 8/9. Moreover, in some cases, it 
may be greater than unity. We refer to this case as the black hole limit. The maximum values 
of H{as,bs; (3) for various values of (3 are shown in Fig.^. Indeed, the maximum mass-to-size 
ratio can sometimes become larger than the black hole limit. However, the physical exterior 
solution generically does not have an event horizon in scalar-tensor theories, in contrast to 
general relativity, and, therefore, the condition, 2M / R > 1, does not imply the existence of a 
black hole. 

Now suppose that a space rocket approaches a star for which 2M/R > 1 and goes into its 
Schwarzschild radius defined by 2M . A spaceman in the rocket would be resigned to his fate 
to die, but we know that he still has a chance to return alive from a false black hole. 



2 \ a. 



(3.47) 



4 Numerical results 



Equations ( |3.6| ) ^ ( |3.9| ) are numerically solved to obtain an interior solution. This solution 
is then matched to the exterior one, and numerical values of the parameters a^, bs and Cs are 
calculated. Some details of the numerical methods are summarized in Appendix C. Since we 
take ifo = 0, G* is equal to G (Appendix A). Hereafter, we use units in which G* = G = 1. 



As for the matter, we assume the following polytropic equation of state |13 

p = Kmtno ( — ] , (4.2) 

nib = 1-66 X 10^24g^ (4.3) 
no = O.lfm"^. (4.4) 

We take the parameters values, F = 2.34 and K — 0.0195 which fit a realistic equation of 
state of high density nuclear matter, and probably also that of a neutron star quite well. Our 
numerical solutions are therefore parametrized by /3 and nc = n(0). It has been shown numer- 
ically that significant effects of ip appear when f3 < —4.35 1^,0, and we are mostly interested 
in cases of negative values of (3. In cases of positive values of we cannot numerically find 
any significantly different behavior of the solutions compared with those in general relativity, 
and any further discussion in these cases is no longer done. 

First, we examine whether our assumption, p'jf < 0, is satisfied. In Fig.Q, we give an 
example of numerical behavior of the effective density for /? = —5 and nc/no = 10. Including 
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this case, we find that the assumption, p[„{r) < 0, is mostly satisfied, as summarized in the 



3rd column of Table 1. Differentiating (3.12), we obtain 



= AH^)p' + AA\^)a{^yp + ( J . (4.5) 



In Fig.|^, we show the 1st, 2nd and 3rd terms in (4.5) with the same parameters as Fig.^. The 



first term is indeed dominant and always negative. In Fig.|6|, we compare the corresponding 
physical quantities, p{r) and 4'{r)^ and the unphysical scalar field, <^(r). It is found that the 
local gravitational constant^ G{4>) = l/4>, increases as p decreases toward the surface. However, 
this behavior is strongly dependent on the coupling function, and the sign of [3 is crucial in the 
present case. 

Next, we examine an extreme example in which the assumption, p'^^ < 0, is violated. Figj^ 
shows the effective density, p^ff, for /? = — 30 and nc/no = 10. It is seen that 

Pcff remains 

constant in the central part and then increases between the two rectangles in FigJ^. In Fig.||, 



we show the three terms in (4^) and find that the positive second term becomes partially 
dominant. Again, this behavior is strongly dependent on the coupling function. In the present 
case, we have a(Lp) = J3(p, where /3 < 0. Therefore, when ip' < and is large, such that 
the second term in ( fl.5| ) is dominant, p'jj becomes positive, and the assumption is violated. In 
Fig.^, we compare p{r), (f>{r) and (p{r). It is found that, despite small values of ip, cj) can be large 
due to a large value of |/3|. However, it should be noted that this assumption is concerned with 
the unphysical quantity, Poff, and that its violation does not necessarily mean that this extreme 
case is unreal. In Fig.^ we show the energy density, p(r), and the pressure, p{r), in the 
physical frame in the extreme case: f3 = —30, Uc/no = 10. The behavior of these quantities 
seems ordinary, that is, they are monotonously decreasing functions of r. Accordingly, one 
may think that this can be a physically acceptable equilibrium solution despite the violation 
of the assumption, p'jj < 0. However, we are forbidden to take (3 smaller than —5 because of 
experimental constraints (Appendix A, 0,(l^). 

For each value of (3, the mass-to-size ratio, 2M/R = H{as, 6^; /3), can be numerically calcu- 
lated as a function of Uc- By changing ric, we search for a maximum value of H{as, bg] (3) for 
each f3. In the 5th column of Table|l|, we summarize our results for the maximum mass-to-size 
ratio, where the parameters are chosen such that the assumption, p'^^ < 0, is satisfied. For 
(3 < —12.07, we find numerically that the assumption is always violated. The first interesting 
example is found in the case, (3 — —12.07, in which the maximal mass-to-size ratio is obtained 
as H-tjiAx = 1.018 when ric/no = 11.2. This is a case in which H^max exceeds the black hole limit, 
H — I. Another interesting example is found in the case, /3 = —11, in which -//max = 0.919 
when nc/riQ ~ 11.3. This is a case in which H^ak exceeds Buchdahl's limit, H — 8/9 w 0.889. 
These examples have academic importance in the sense that our analytic results in the previous 
section are partially realized also in the numerical solutions. 

To this point, the stability of our numerical solutions has not been taken into account. We 
have found that p — Zp may be a good estimator of the stability, as described below. The 
baryonic mass of a star is defined as |^ 

/ 2m \ ~^ 

m^mb AT:nA^{ip)r^ 1 dr. (4.6) 







We have numerically examined how ffi depends on Uc and find a significant correlation between 
the signature of drh/dric and that of p — 3p. That is, the cases, p — 3p < and p — 3p > 0, 
approximately correspond to the cases, dm/dric < and dfh/dric > 0, respectively. Accord- 
ingly, we shall interpret the signature of p — 3p as a measure of the onset of the instability in 
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our numerical calculations. When we impose the condition, p — 3p > 0, we cannot numerically 
find cases in which Hmax exceeds Buchdahl's limit, as is seen in the 4th and 5th columns in 
Table |l]. Note that the condition, /3 < —5, also excludes all the interesting cases in which H^ax 
exceeds Buchdahl's limit. 

Let us find a critical value, /3c < 0, of (3, such that, for f3 < /3c, nonlinear behavior of the 
scalar field begins to appear. We numerically calculated a^, &s and Cg as functions of ric for 
/3 = —4, —5 and —6 under the conditions, p'^^ < and p — Sp > 0. We show (a^, bg) and Cs 
in Fig.|l^ and Fig, 12, respectively. For /3 > —4, almost no deviation from general relativity 



appears. In the cases that (3 = —5 and —6, these parameters show deviations from general 
relativity in which as = bg and Cg — 0. Our results are consistent with the previous works [^0, 
in which /3c is found to be —4.35. Note that, even when /3 < /3c, our inequality, \cs\ < 2v3/9, 
is surely satisfied. This reconfirms our assertion that the nonlinear effects are always bounded 
in this sense. 

Now we shall briefly compare our numerical results with those in previous works 1^ , . The 
maximum baryonic mass of a star is defined as the peak of the fh-ric relation. We numerically 
found that the maximum baryonic mass increases from the general relativistic value, 2.23Mq, 
to 2.38M0 and 2.96M0 for /3 = —5 and —6, respectively. The corresponding radius defined 



by (3.43) also increases from the general relativistic value, 11.0 km, to 12.0 km and 12.9 km 
for /3 = — 5 and —6, respectively. In Ref.|^, a fractional binding energy, /be = 2m/6— 1, is 
used as a measure of the scalar field contribution to the mass. We numerically found that the 
maximum value of /be increases from the general relativistic value, /be = 0.14, to /be = 0.16 
and 0.22 for f] — —5 and —6, respectively. Though a shghtly different asymptotic value of (po 
has been adopted in Ref.||^, these results are consistent with the previous results. 

The mass-size relation of neutron stars has been thoroughly studied in general relativity 
by solving the Oppenheimer-Volkoff equation, and it has been found numerically that, as 
the equation of state becomes softer, the mass-to-size ratio becomes larger when its mass is 
fixed|l|. In Fig.^ we compare, under the condition, p — 3p > 0, the relations between the 
mass-to-size ratio and the mass in general relativity and in scalar-tensor theories. It is seen in 
Fig.|l^ that the deviation from general relativity due to the scalar field begins appearing for 
M > 1.2M0. When M < I.TMq (M > l.TM©), our numerical solutions in the scalar-tensor 
theory correspond to the solutions in general relativity with the softer (stiffer) equation of state. 
The mass of PSR1913-I-16 has been evaluated as 1.4M0 [13], and, if the adopted equation of 
state is adequate, the scalar field contribution to the mass-to-size ratio is negligibly small when 
/3 > — 5. Further discussion on the equation of state is left as a future work. 

Finally, we shall derive a redshift formula in the scalar-tensor theories. A null vector, k'^, 
tangent to the radial null geodesic, and a four- velocity, C/^, of a static observer are, respectively, 
given by 

^ A-^{e-'^, 1,0,0), C/^ = (A-^e-^/^ 0,0,0), C/^C/^ = -1. (4.7) 
The frequency, w, of a light ray is given by 

CO ^ -g^.k'^U'' ^ A-^e"^/'. (4.8) 
The redshift, z, is then obtained as 



(4.9) 



where the observer is assumed to be at the spatial infinity, a/x — > 0. By using (3. 22), (3. 23) 
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and (3.40), we obtain the redshift formula in the present specific scalar-tensor theory as 



l + z = (l-a,)-^expi--/3 ^ln(l-a,) 



2 



1 




(¥) 



The maximum value of z depends on /3 and the parameters, a a and hs, in the allowed region, 
D, in Fig.l. Theoretically, the possible maximum value of Zmax, is obtained as z^^x = 2 
in general relativity and z^^^ = 164 and 356 for /3 = — 5 and (3 = —6, respectively. We 
compare these values with those obtained in the numerical calculations under the condition, 
p — Bp > 0. We numerically find that z^^x = 0.43 in general relativity, and that z^^x = 0.44 and 
0.55 for /3 = — 5 and —6, respectively. In our numerical solutions with the presently adopted 
model parameters, the deviation of the redshift formula from that in general relativity remains 
small compared with the theoretically possible deviation. However, the redshift difference, 
0.01, numerically found for f3 = —5 may be detectable. Therefore the redshift measurement of 
neutron stars will provide us with a possible tool for the experimental test of general relativity. 



5 Summary 

We have derived a modified Buchdahl inequality in scalar-tensor theories of gravity. As 
a result, we have obtained two theory-independent inequalities, bg < 8/9 and \cs\ < 2\/3/9. 
The first inequality corresponds to the Buchdahl inequality in general relativity. The second 
inequality is characteristic of scalar-tensor theories. Consequently, even if the scalar field is 
locally amplified due to non-perturbative effects in a strong gravitational field, the characteristic 
amplitude of the scalar field, \cs\, is bounded in this sense. 

The modified Buchdahl inequality is then reformulated to obtain a theory-dependent mass- 
to-size ratio, 2M/R, with an example of the coupling function, A{(p), in a simple form. If we 
take tpo = 0, the mass in the physical frame is the same as that in general relativity, M = 6/2. 
However, the physical radius, R, of the star can be smaller than the general relativistic one. 
As a result, the mass-to-size ratio can exceed not only Buchdahl's limit but also the black hole 
limit in contrast to general relativity. 

Our analytic results have been numerically confirmed when we assume a polytropic equation 
of state for the matter. In particular, we have found numerical solutions in which the mass-to- 
size ratio exceeds both Buchdahl's limit and the black hole limit. However, these theoretically 
interesting stars could not be found numerically under the condition, p — 3p > 0, which is 
interpreted as a numerical measure of the onset of the instability of a star. Moreover, under 
this condition, we find numerically that any quantitative deviation from general relativity 
due to the scalar field remains comparatively small in contrast to our analytic results, where 
possible significant effects of the scalar field are expected. However, as discussed briefly, some 
measurable effects in astronomical observations may exist. 

Now suppose that a space rocket approaches a massive star for which 2M/R 1. If the 
rocket accidentally goes into a Schwarzschild radius of the star defined by 2M, a spaceman 
in the rocket would be resigned to his fate to die. Now we know that, unfortunately for him, 
even if scalar-tensor theories describe classical gravity, he would hardly have a chance to return 
alive, because he could hardly meet a real false black hole. He has two possible futures, and 
they are equally tragic: 

• If it is a black hole in terms of general relativity, he can never escape. 
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• If it is a naked singularity in terms of scalar-tensor theories, nobody knows what will 
happen when he touches it. 
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A Observational constraints 



In general, the coupling strength, a{ip), can be an arbitrary function of ip, and in the limit, 
a((p) 0, scalar-tensor theories approach general relativity. One defines 

ao = a{ipo), (A.l) 



dip 



(A.2) 



where ipo is the asymptotic value of ip at spatial infinity. In the post-Newtonian approximation, 
the PPN parameters and the effective gravitational constant are expressed as follows [0 : 

^-lE = (A.3) 
G = G,A2((^o)(l + a2((^o)). (A.5) 



Pe-1 - :,7f^, (A.4) 



General relativity corresponds to the case that /3b = 7_e = 1 @l,(l5|. Experiments on the time 
delay and deflection of light in the solar system constrain |1 — 7£;| as Q 

|1 - 7^1 <2x 10-3, (A.6) 

which constrains ijj{4>) and as 

u > 500, al < 10-3. (A.7) 

The lunar-laser-ranging experiments constrain j/Je — 1| as JT^ 

l/Sfi-ll <6x 10-^ (A.8) 

which only constrains some combination of ao and Po- Consequently, if ao tends to zero, 
the constraint on (3q is effectively lost. However, by adopting a specific coupling function, 
A = exp (i/3(^^), another constraint on /3o is obtained from observations of the binary-pulsars, 
PSR1913-hl6, as §,|l4| 

f3o > -5. (A.9) 



When we take (3.4C) as a coupling function, the coupling strength is a{ip) — f3ip. Ac- 
cordingly, ao — f3ipo, and we obtain the constraint on ipo. In this paper we take ipo ^ for 
simplicity. 
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B An exterior solution 



A line element of the Einstein frame in the Just coordinate is ||l6| 

ds^ = -e<^^Ut^ + e-T(x)dx^ + e^^^^-''^^^dn'^. (B.l) 
The field equations in the exterior space-time are 

y + YA' = 0, (B.2) 

-7'' + 7'A' - + 7" - 2A" = V^, (B.3) 

2 + e^(7'A'-A'2 + 7"-A") = 0, (B.4) 

(p" + AV' = 0, (B.5) 



where a prime denotes differentiation with respect to %. With ( |B.2| ),(B.4) and (B_^), the 
exterior solution can be obtained as 

e^(x) = x'-ax, (B.6) 

(^(X) = (po + -lnfl--V (B.8) 

where a, b and c are constants of integration, and (^q denotes the asymptotic value of Lp at 
infinity. With (B.3), one finds 

- 52 = 4c2. (B.9) 

The coordinate transformation between the Schwarzschild coordinate, r, and the Just co- 
ordinate, X, is given by 

a-b 

Note that r ^ x &t spatial infinity. In the Schwarzschild coordinate, a line element becomes 

dsl = -e^^'Mdt^ + e^^Mdr^ + r'^dn'. (B.ll) 
The exterior solution in the Schwarzschild coordinate is given by 

Asymptotic behavior of the exterior solution at spatial infinity are as follows: 

e^'^^'^) — ^ 1--, (B.14) 
r 

e2^('-) ^ 1 + (B.15) 
r 

<p(r) ^ (^0--. (B.16) 
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C An interior solution: Numerical methods 



Using the variables, TO»(r) and defined by 



the field equations 
dm* 



^ (^.S|) become 



K{r) 



2m*(r) 



1 



dr 
dv 
dr 
dip 
dr 

dr 
dp 
dr 



, ^ +lry = <^{r), 
r(r — 2,m^) I 



+ 47rG'*A^(i^)r^p 1 



a(^)(p - 3p) + (p-p)rV'] 



r — 2m, 
-(p + p) ($ + a((p)V') . 
The total baryon mass measured in the physical frame is 

,0j3 



m = mb / n^—gu d x — nif, / ATrnA {ip)r 1 — 



2(r — m*) 
r(r — 2m*) 



V', 



2m, 



dr. 



(C.l) 

(C.2) 
(C.3) 
(C.4) 
(C.5) 
(C.6) 

(C.7) 



Given the equation of state, we can numerically integrate the above field equations outward 
from the center, r — 0, with the boundary conditions as follows: 



(C.8) 



m»(0) 


- 0, - 


(^(0) 


= Vc, 


m 


- 0, > 


p(0) 


= Pc, 


P{0) 


= Pc, , 



where pc and pc are given by replacing n in (41) and ( [4. 4) with rtc = n(0). The surface of 
a star, r = r^, is determined by the condition, pirg) = 0. A numerically obtained interior 
solution is to be matched to the exterior one by the conditions |^ , (t) : 



: tanh 



b = 



2r 



2m, 



y v's + l/rs 



exp 



I'. 



: tanh 



l/rs 



2v' 



Vfo^ + 4c2, 



(C.9) 

(C.IO) 

(C.ll) 
(C.12) 



where a prime denotes differentiation with respect to r, and the subscript, s, refers to quantities 
evaluated at the surface, Ts. The central value of Lp, ip^^ is chosen such that we have Lpo = 0. 
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Figure 1: The allowed region, 13, is depicted, where horizontal and vertical axes denote, re- 
spectively, bs and as. The characteristic points, P = (4(3 - 2^2), 4(3V2 - 4)), Q = (8/9,8/9) 
and K = (4/9,8/9), are shown. In general relativity, — bg. Buchdahl's limit is denoted by 
Q. On K, \cs\ takes the maximum value, |cs|max = 2-\/3/9. 



Figure 2: On each hues, H{as,bs; (i) — 8/9 for various values of /3: 0, —3, —5,-10 and —100. 
Horizontal and vertical axes denote, respectively, bg and a^. 



Figure 3: The Maximum mass-to-size ratio, Hmax, is shown as a function of 13. Horizontal 
and vertical axes denote, respectively, /3, and TJmax- A horizontal line, Hj^j^^ — 8/9, denotes 
Buchdahl's limit in general relativity. When (3 < 0.4, Hmax exceeds Buchdahl's limit. When 
/3 < 0.2, i?MAx exceeds unity, a black hole limit. 



Figure 4: The effective density, p^{{{t), is shown in the case that (3 — —b and nc/no — 10. 
Horizontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, 
and the effective density, p^ff(r)/(mf,no). It is seen that the assumption, jo'jj(r) < 0, is satisfied. 



Figure 5: Each term of ( |4.5| ) in p'^^fir) is shown for (3 = —h and nc/riQ = 10. Horizontal and 
vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, and the 1st, 
2nd and 3rd terms in ([1.5|). 



Figure 6: We compare 0(r), p{r) and Lp{r) in the case that (3 = —5 and nc/riQ = 7.9. Hori- 
zontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, and 
p{r) / (uibno) , (j){r) and ip{r). On a thin dotted line, G = l/4>() — 1, i.e., general relativity. 



Figure 7: The effective density, p^if{r), is shown in the case that (3 = —30 and ric/no — 10. 
Horizontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, 
and the effective density, p^!f{r) / {rribno) . It is seen that the assumption, p[ff{r) < 0, is partially 
violated between two rectangles. 



Figure 8: Each term of ( |4.5| ) in p[„{r) is shown for f3 — —30 and ric/no = 10. Horizontal and 
vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, and the 1st, 
2nd and 3rd terms in ([l.5|). 



17 



Figure 9: Wc compare (/>(r), p{r) and (p{r) in the case that /? = —30 and Uc/uq = 10. Hori- 
zontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, and 
p(r)/(m(,no), (p{r) and <^(r). On a thin dotted line, G = l/(po = 1, i-e., general relativity. 



Figure 10: We show p{r) and p{r) in the physical frame for (3 = —30 and Uc/no = 10. 
Horizontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, 
and p{r) / {mbTio) and p{r) / {rribno) . Though the assumption, P3ff(r) < 0, is violated in the 
Einstein frame, the conditions, p(r) — 3p(r) > 0, p'{r) < and p'{r) < are all satisfied in the 
physical frame. 



Figure 11: We show the parameters, (0^,65), in each equilibrium solution for ric/rio = 2.5 ~ 
10.3. We take (3 = —4, —5 and —6, and impose the conditions, p{r) — 3p(r) > and p^jf (r) < 0. 
Horizontal and vertical axes denote, respectively, bg and a^. 



Figure 12: We show the parameter, c^, in each equilibrium solution for ric/no = 1.0 ~ 10.3. 
We take = —4, —5 and —6, and impose the conditions, p{r) — 3p{r) > and p^ff(r) < 0. 
Horizontal and vertical axes denote, respectively, ric/no and Cg- A horizontal thin dotted line 
denotes the limit on Cg, i.e., Cg = —2y/3/9. 



Figure 13: We show the relation between the mass-to-size ratio and the mass. We take /3 = — 5 
and —6, and impose the conditions, p{r) — 3p{r) > and plff{r) < 0. Horizontal and vertical 
axes denote, respectively, M/Mq, and 2M/R. The solid line represents the mass-to-size ratio 
in general relativity. 
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/3 


ric/no 


p:« < 


p-3p>0 


H 


-12.07 


0.1 ~ 10.3 








0.679 




10.3 - 11.2 





X 


1.018 


-11.0 


0.1 ~ 10.3 








0.666 




10.3 ~ 11.3 





X 


0.919 


-10.0 


0.1 ~ 10.3 








0.651 




10.3 ~ 11.4 





X 


0.834 


-6.0 


0.1 ~ 10.3 








0.556 




10.3 - 11.2 





X 


0.569 


-5.0 


0.1 ~ 10.3 








0.517 




10.3 ~ 15.7 





X 


0.567 


-4.0 


0.1 - 10.3 








0.514 




10.3 ~ 14.9 





X 


0.562 


0.0 


0.1 ~ 10.3 








0.514 




10.3 ~ 14.9 





X 


0.562 



Table 1: We summarize our numerical results. In the 1st column, /3 is given. In the 2nd 
column, we give a range of ric from our numerical studies. In the 3rd and 4th columns, we 
indicate, respectively, whether the assumption, p^^j < 0, and the condition, p — 3p > 0, are 
satisfied. In the 5th column, the maximum mass-to-size ratio is shown for each (3. 
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Figure 1: The allowed region, D, is depicted, where horizontal and vertical axes denote, re- 
spectively, bs and a^. The characteristic points, P = (4(3 - 2i/2), 4(3\/2 - 4)), Q = (8/9,8/9) 
and K = (4/9,8/9), are shown. In general relativity, = bg. Buchdahl's limit is denoted by 
Q. On K, \cs \ takes the maximum value, |cs|max = 2-\/3/9. 
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Figure 2: On each lines, H{as, bg] (3) = 8/9 for various values of (3: 0, —3, —5, —10 and —100. 
Horizontal and vertical axes denote, respectively, bg and a^. 



1 




Scalar tensor theory 



General relativity 



1.5 



2.5 



Figure 3: The Maximum mass-to-size ratio, ff, 
and vertical axes denote, respectively, /3, and H, 
Buchdahl's limit in general relativity. When < 0.4 
/3 < 0.2, i?MAx exceeds unity, a black hole limit. 



is shown as a function of /3. Horizontal 
MAX- A horizontal line, H^max = 8/9, denotes 



exceeds Buchdahl's limit. When 




0.6 

Radial coordinate 



Figure 4: The effective density, Pcff{r), is shown in the case that /3 = —5 and ric/no = 10. 
Horizontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, 
and the effective density, p^Hif')/{'mbno)- It is seen that the assumption, p'ff(j') < 0, is satisfied. 
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Figure 5: Each term of (4.5) in /9'„(r) is shown for (3 = 
vertical axes denote, respectively, the radial coordinate, 
2nd and 3rd terms in (4.5). 
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Figure 6: We compare 4>{r), p(r) and ip{r) in the case that (3 = —5 and tIc/uq = 7.9. Hori- 
zontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, and 
p(r)/(mf,no), 4>{r) and ip{r). On a thin dotted line, G = l/4>o = 1, i-e., general relativity. 
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Figure 7: The effective density, p^if), is shown in the case that /3 = —30 and ric/no = 10. 
Horizontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, 
and the effective density, Ppff(r)/(mf,no). It is seen that the assumption, p'ff(j') < 0, is partially 
violated between two rectangles. 
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Figure 8: Each term of (4.5) in p'(j(r) is shown for /3 = —30 and ric/no = 10. Horizontal and 
vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, and the 1st, 
2nd and 3rd terms in (4.5). 
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Figure 9: Wc compare (j){r), p{r) and ^p{r) in the case that [i = —30 and ric/no = 10. Hori- 
zontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, and 
p(r)/(mf,no), 4>{r) and ip{r). On a thin dotted line, G = l/4>o = 1, i-e., general relativity. 
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Figure 10: We show p(r) and p{r) in the physical frame for f3 = —30 and Hc/hq = 10. 
Horizontal and vertical axes denote, respectively, the radial coordinate, r, in the unit of 10km, 

and p{r)/(m,[,n{)) and p{r) / {mf,no) . Though the assumption, p[.ff{r) < 0, is violated in the 
Einstein frame, the conditions, p{r) — 3p{r) > 0, p'(r) < and p'{r) < are all satisfied in the 
physical frame. 
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Figure 11: We show the parameters, {as,bs), in each equilibrium solution for ric/no = 2.5 ~ 
10.3. We take (3 = —4, —5 and —6, and impose the conditions, p(r) — 3p(r) > and p'ff(j') < 0. 
Horizontal and vertical axes denote, respectively, and a^. 
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Figure 12: We show the parameter, c^, in each equilibrium solution for n^/no = 1.0 ~ 10.3. 
We take /3 = —4, —5 and —6, and impose the conditions, p{r) — 3p(r) > and p'ff(j') < 0. 
Horizontal and vertical axes denote, respectively, ric/no and c^. A horizontal thin dotted line 
denotes the limit on c^, i.e., = — 2-\/3/9. 
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Figure 13: We show the relation between the mass-to-size ratio and the mass. We take /3 = —5 
and —6, and impose the conditions, p(r) — 3j)(r) > and p^„(^) < 0. Horizontal and vertical 
axes denote, respectively, M/M.Q, and 2M/R. The solid line represents the mass-to-size ratio 
in general relativity. 
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Table 1: We summarize our numerical results. In the 1st column, f3 is given. In the 2nd 
column, we give a range of ric from our numerical studies. In the 3rd and 4th columns, we 
indicate, respectively, whether the assumption, p'jj < 0, and the condition, p — 3p > 0, are 
satisfied. In the 5th column, the maximum mass-to-size ratio is shown for each f3. 
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